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1 .  INTRODUCTION 

Ships  and  submarines  frequently  radiate  sounds  underwater  which  exhibit 
distinct  amplitude  modulation.  For  example,  the  noise  from  a  cavitating 
propeller  is  usually  modulated  at  a  frequency  equal  to  the  blade  rate  (i.e., 
the  rate  of  rotation  of  the  shaft  multiplied  by  the  number  of  blades  of  the 
propeller).  This  noise  can  be  received  by  a  passive  sonar,  and  with  training 
and  in  conjunction  with  other  sounds,  an  operator  can  often,  from  the  modulation, 
identify  the  type  of  vessel  and  estimate  its  speed. 

The  operator  can  be  greatly  assisted  by  rectifying  the  sonar  signals  to 
extract  the  modulation  and  spectrum  analysing  the  result.  This  is  termed 
"DEMON"  processing.  A  block  diagram  of  a  DEMON  processor  is  shown  in 
figure  1 .  The  writer  is  not  aware  of  any  published  mathematical  analysis 
of  this  processor.  In  this  memorandum,  an  expression  of  its  theoretical 
performance  is  derived. 

It  should  be  pointed  out  that  the  processor  shown  in  figure  1  is  not 
necessarily  the  optimum  detector  for  such  modulated  signals.  Tuteur(l)  has 
made  an  attempt  to  determine  the  optimum  processor  and  its  performance,  but 
the  simplifying  assumptions  he  makes  are  unrealistic,  and  as  a  consequence  his 
results  should  be  regarded  as  useful  only  in  giving  an  upper  bound  to  the 
performance  which  might  be  achieved.  The  present  writer  has  not  yet  succeeded 
in  determining  the  optimum  processor  under  realistic  assumptions. 


2.  ASSUMPTIONS 

The  following  assumptions  are  made. 

(i)  The  sound  radiated  from  the  vessel  has  the  properties  of  zero-mean 
G-aussian  noise  which  is  amplitude  modulated  by  a  low-frequency 
sinusoidal  wave.  The  resultant  is  henceforth  referred  to  as  the 
"signal". 

(ii)  The  frequency  of  the  modulating  wave  is  known,  but  its  phase  is 
random. 

(iii)  The  signal  is  contaminated  by  additive  zero-mean  G-aussian  noise. 

(iv)  Both  the  signal  and  noise  are  wholly  contained  in  a  narrow  bandwidth 
B  (i„e.,  the  centre  frequency  is  much  greater  than  B) . 

(v)  The  spectrum  is  computed  for  samples  of  data  of  duration  T,  where  BT 
is  an  integer  >>  1 . 

(vi)  The  average  is  taken  of  M  spectra,  where  M  »  1. 

(vii)  Signal  and  noise  are  stationary,  ergodic  processes. 

(viii)  Signal  and  noise  are  independent. 

Although  it  is  unlikely  that  the  modulating  frequency  would  be  known  a 
priori,  assumption  (ii)  above  is  nonetheless  not  unrealistic,  because  in 
practice  one  would  compute  the  entire  power  spectrum  over  all  the  frequencies 
of  interest. 
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3.  NOTATION 

Let  the  input  waveform  to  the  processor  be  described  by 

2nPt 


z(t)  =  1  +  b  cos  f  T  +  <pj  s(t)  +  n(t)  , 


(i) 


where  s(t)  and  n(t)  are  zero-mean  Gaussian  random  waveforms, 

P  is  an  integer  >  0, 

P  <<  BT  (assumption  (i)  ), 

0  <  b  ^  1 , 

9  is  random  in  (0,  2tt)  . 

It  is  possible  to  express  s(t)  and  n(t)  in  terms  of  low-frequency  (complex) 


waveforms  modulating  a 

"carrier1 

s(t)  = 

i  i 

V2 ; 

r«(t) 

n(t)  = 

i  i 

V2 ; 

^(t) 

(2) 


(3) 


where  a(t)  and  /3(t)  are  in  general  complex,  and 
*  denotes  the  complex  conjugate. 

We  shall  decompose  a(t)  and  /3(t)  into  their  respective  Fourier  components: 

BT 


a(t)  o:  x(n)  exp(i2w  nt/T) 

n=1 


(4) 


BT 


£(t)  n:  ^  y(n)  exp(i2?r  nt/T)  , 


n=1 


1 

where  x(n)  =  a(t)  exp(-i27r  nt/T)  dt, 
i 


y(n)  =  ~  f  p( t)  exp(-i2 n  nt/T)  dt. 


(5) 


(6) 
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For  BT  »  1, 


<x(n)  x*(m)>  S(n)  8^  , 

<y(m)  y*(®)>  -  N(n)  5^  , 


(7) 


where  S(n)  is  the  unmodulated  signal  power  component  at  the  frequency  corres¬ 
ponding  to  the  index  n, 

and  5  is  the  Kronecker  delta, 

nm  . 


=  0  ,  n  +  m, 

=  1  ,  n  =  m. 


(8) 


N(n)  is  the  noise  power  component  at  the  frequency  corresponding  to 
the  index  n, 

^  )>  denotes  an  ensemble  average. 

By  virtue  of  assumption  (viii), 

<^x(n)  y*(m)^)>  =  0,  for  all  n,  m, 

and  from  assumption  (iv) 

x(n)  =  y(n)  =  0  for  n<  1  and  n  >  BT 


(9) 

(10) 


In  what  follows,  we  shall  make  extensive  use  of  the  following  theorem  for  the 
moment  of  complex  gaussian  processes (2)  : 


<'z(mi)  z(m2)...z(m  )  z*(nO  z*(n2) . .  .z*(n  o)'>  =  0  if  r  *  s, 

^  r  s 

<zK(i))  z*(nt)>  <zK(2))  z*(“2)> <zK(s))^(ns)> 

if  r  =  s  (11) 

where  n(  )  denotes  a  permutation. 

Specifically,  if  r  =  s  =  2, 

<z(m1)  z(m2)  z*(ni)  z*(n2)>  =  <z(m1)  z* (n1 )> <z (m2)  z*(n2)> 

+  <z(m1)  z*(n2)><z(m2)  z*(ni)>  (12) 


4.  MATHEMATICAL  DEVELOPMENT 

After  squaring  and  low-pass  filtering  z(t),  .only  terms  remaining  are  those 
not  containing  f  ,  and  we  have,  from  equations  (1)  to  (3): 
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+  b  cos  (^Tp-  +  95)  J  |a(t)  |2  +  |/?(t)  |2 
+  jj1  +  b  cos  (^Tp-  +  4)  j  |4(t)/?*(t)  +  a*(t)/9(t)j 


After  substituting  equations  (A)  and  (5)  into  this  expression,  and  a  little 
manipulation,  we  have: 


b2  ^  /i27rPt  .  /  i2mPt 

1  +  -r  +  t  exp  — - —  +  i<p  +  b  exp - - 1  9 


+  exp 


iAflPt 


iAflPt 


+  129  +T-  exp  -  i2<p  . 


^4  x(n)  x*(m)  exp  |^i27r(n-m)  t/T 


n,m 


y(n)  y*  (m)  exp  jjL27r(n- 


m 


)  t/T 


J 


n,m 


,  „  b  /i27rPt  .  \  b  /  i27rPt  .  \ 

+  1 1  +  ~2  exP  r — r —  +  19 J  "2  exp  ( - - t - 19y  j 


n,m 


x(n)  y*(m)  exp  j  i27r(n- m)t/T  j  +  x*(n)y(m)  exp  \  -i27r(n-m)  t/T 

(13) 


Consider  now  multiplying  this  by  cos(27rPt/T)  and  integrating  over  (0,T). 
Noting  that 


Y  f  exp  fi27r(n- m)  t/T*)  dt  =  , 


nm 


(14) 
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we  have 


i  1  + 


x(n)  x*(n  +  P)  +  x(n)  x*(n-P) 


be 


1<P 


x(n)  x*(n+2P)  +  x(n)  x*(n) 


190  v— » 


+  \  |  x(n)  x*(n)  +  x(n)  x*(n-2P) 


+ -vr-  e'^2'’0  'V  fx(n)  x*(n+3P)  +  x(n)  x*(n  +  P) 


-jr  e  i2<P  fx(n)  x*(n-  P)  +  x(n)  x*(n-  3P) 


+  ^  |y(n)  y*  (n  +  P)  +  y(n)  y* (n  -  P) 


+  %  ^|*(n) 


n+P)  +  x(n)  y*(n-P) 


+ 


_b 

4 


+ 


b. 

4 


+ 


JL 

2 


+ 


b_ 

4 


^  X!  [x(n)  y*^n+2p^  +  x(n)  y*(n)j 


e  ^  |x(n)  y*(n)  +  x(n)  yii<(n-  2P) 


x*  (n)  y(n-P)  +  x*(n)  y(n  +  P) 


} 


e±(p  \  x*  (n)  y(n  -  2P)  +  x*(n)  y(n) 


e  1<P  x*(n)  y(n)  +  x*(n)  y(n  +  2P) 


(15) 


A  similar  expression  (except  for  a  change  of  sign  inside  each  curly  bracket, 
and  the  presence  of  i  outside  each  term)  obtains  after  multiplying  by 
sin(27rPt/T)  and  integrating  over  (0,T). 
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These  two  outputs  (i.e.,  after  cosine  and  sine  multiplications  and  integrating) 
are  squared  and  summed.  The  algebraic  manipulations  are  simplified  if  we  note 
that  the  two  outputs  take  the  forms: 

cos :  (p  +  q)  +  (r  +  s)  +  . . . 

sin:  (ip  -  iq)  +  (ir  -  is)  +  ..., 

and  that  after  squaring  and  summing  the  only  terms  remaining  are  the  cross 
products : 

4  pq  +  4  ps  +  4  qr  +  4  rs  +  . . . 

Making  use  of  equations  (7)5  (9)  and  ( 1 2 ) ,  we  find  that  the  expression  for 
the  mean  after  squaring  and  summing  eventually  reduces  to: 


1  +  tT  y  S(n)S(n  +  P)  +  Y  S(n)s(n+  3P)  +  *£  )  S(n)S(n-  P) 


+  b2  \  S(n)S(n  +  P)  +  b2  >  S2(n)  +  b2  \  S(n)  +  \  N(n)N(n  +  P) 


+  \  S(n)N(n  +  P)  + 


S(n)N(n-P)  +  ■— 
'4 


S(n)N(n+  2P)  +  )  S(n)N(n-  2P) 


+  y-  )  S(n)N(n) 


(16) 


In  expression  (16),  the  sums  are  taken  over  the  range  (1,BT);  it  should  be 
noted  that  from  the  initial  assumption  (iv), 

S(n)  =  N (n)  =  0  for  all  n  <  1  and  n  >  BT. 


We  now  make  the  additional  assumption  that  except  at  the  band  limits,  both  the 
signal  and  noise  spectra  are  locally  slowly  varying,  i.e., 

S(n  +  2P)  -  S(n  +  P)  *  S(n), 

and  N(n  +  2P)  -  N(n  +  P)  a  N(n). 

Expression  (16)  for  the  mean  output  is  then  approximated  by 


b2 


2 


+ 


+  (2  +  b2) 


(n)N(n)l 


(17) 
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We  are  also  interested  in  the  variance  of  the  output  when  signal  is  absent. 
In  the  absence  of  signal,  S(n)  =  0,  and  the  mean  output  is  simply 


The  output  without  signal  is 


^  y(n)  y* (n  +  P)  y(m)  y*(m  -  P) 
n,m 


The  mean  square  output  is 


<y(k)  y*(k  +  P)  y(0  y*(«-P)  y(m)  y*(m  +  P)  y(n)  y*(n-P)>  (18) 

k,  ^m,n 

Using  equations  (7)  and  (ll),  we  find,  after  considerable  algebraic  manipulation, 
that  the  mean  square  output  for  noise  only  is 

2  +  6Ij[N(n)]4  (19) 

The  variance  of  the  output  is  therefore 

mean  square  output  -  (mean  output)2  = 


The  signal-to-noise  power  ratio  at  the  output  of  the  processor  may  be  defined 
to  be: 


N(n) 


+  6  \  N(n)  (20) 


d 


£  (change  in  mean  with  signal  present)2 
~  variance  without  signal  — 


As  we  add  M  independent  spectra,  it  can  readily  be  seen  that  the  signal-to-noise 
output  power  ratio  will  be 
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For  the  special  case  in  which  both  (unmodulated)  signal  and  noise  have  flat 
spectra,  let 


Defining 


S(n)  =  S,  N(n)  =  N  for  all  n 


K  =  BT, 


N  jJ 


KM  fb2  K  ju2  +  (l  +  3b2  +-g-b4")  jj.2  +  (2  +  b2)  n 


K  +  6 


In  most  cases  of  practical  interest,  K  >>  1 ,  (J  «  1 .  Hence 


we  have 


d  -  M^2  b2  Kju  +  (2  +  b2) 


If  M  is  sufficiently  large,  the  output  will,  by  the  central  limit  theorem,  be 
approximately  gaussian,  and  could  therefore  be  defined  by  second-order  statistics. 
Further,  we  could  then  use  curves  for  Receiver  Operating  Characteristics  (as 
derived,  for  example,  in  references  3  and  4)  directly;  these  give  a  plot  of 
the  probability  of  detection  versus  the  probability  of  false  alarm  for  various 
values  of  d2  (see  figure  2),  which  in  our  case  is  simply  the  output  signal-to- 
noise  ratio. 

It  should  be  noted  that  in  equations  (21)  to  (23)  ju  is  not  the  input  signal- 
to-noise  ratio,  because  S  is  the  power  spectral  density  of  the  unmodulated 
signal.  The  modulation  will  increase  the  actual  input  signal-to-noise  power 


ratio  to  a  value  of  +  — mJ  (J  . 

It  is  also  noteworthy  that  the  results  obtained  here  are  considerably  more 
pessimistic  than  those  obtained  by  Tuteur(l),  who  assumes  knowledge  of  the 
phase  of  the  modulating  signal.  Essentially,  his  output  signal-to-noise  ratio 
varies  with  ju2,  whereas  the  result  derived  here,  which  assumed  random  phase, 
varies  with  ju4  for  large  K. 

5.  COMPARISON  WITH  SIMPLE  ENERG-Y  DETECTOR 

By  way  of  comparison,  consider  the  simple  energy  detector  shown  izi  figure  3- 
The  output  signal-to-noise  ratio,  defined  as  before  to  be 

^  a  (change  in  mean  with  signal  present)2  g  ^ 

variance  without  signal  ^  ^  \fKo  v v 

may  be  shown,  after  the  usual  algebra,  to  be 


h2\  2 

1  +  id2  BTM 
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Comparing  equations  (23)  and  (24),  we  find  that  for  the  Demon  processor,  with 
small  signals,  the  output  signal-to-noise  power  ratio  varies  as  [iA ,  whereas  with 
the  simple  power  detector  it  varies  as  p2. 

Hence  it  is  to  be  expected  that  in  general  detection  using  a  simple  power 
detector  would  take  place  before  Demon  classification  could  be  effected. 


6.  EXAMPLE 


By  way  of  a  simple  example,  consider  the  case  of  noise  and  (unmodulated) 
signal  spectra  which  are  flat.  It  is  required  to  calculate  the  value  of  p 
which  will  give  a  probability  of  detection  of  90/?,  and  a  probability  of  false 
alarm  of  001/?.  Other  given  parameters  are: 

B  =  1 03  Hz  , 

T  =  10s, 

Mi  =  100, 

b  =  1-0 

From  figure  2,  it  is  seen  that  d  =  25.  Hence 

100  ii2  [  104  (i  +  3j2  =  25 

i.e.  ,  p  ^  7  x  10'3  , 

and  the  actual  signal-to-noise  input  ratio  is 


or  about  -20  dB. 

Under  the  same  circumstances,  the  ordinary  power  detector  of  broad-band 
signals  would  have  a  detection  threshold  of  -25 dB. 


7.  CONCLUSIONS 

An  expression  has  been  obtained  for  the  output  signal-to-noise  ratio  of  the 
Demon  processor.  It  is  found  that  in  general  the  detection  (using  conventional 
power  detection)  would  occur  before  Demon  classification  could  be  effected. 

For  small  signals,  the  signal-to-noise  power  ratio  out  of  a  Demon  processor 
falls  as  the  fourth  power  of  the  input  signal-to-noise  ratio,  whereas  that  of 
the  power  detector  falls  as  its  square.  This  means  that  the  threshold  of  the 
Demon  processor  can  be  expected  to  appear  more  sharply  defined  than  that  of  the 
power  detector. 

As  would  be  expected,  the  performance  of  the  Demon  processor  falls  rapidly 
as  the  depth  of  modulation  of  the  signal  is  reduced. 
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Figure  2.  Receiver  operating  characteristics 
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